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We note that the theoretical prediction of neutrinos from gamma-ray bursts (GRBs) by IceCube 
overestimates the GRB neutrino flux, because they ignore both the energy dependence of the fraction 
of proton energy transferred to charged pions and the radiative energy loss of secondary pions and 
muons when calculating the normalization of the neutrino flux. After correction for these facts the 
GRB neutrino flux is reduced, e.g., by a factor ~ 5 for typical GRB spectral parameter, and may 
be consistent with the present zero event detected by IceCube. More observations are important 
£\J \ to push the sensitivity below the prediction and test whether GRBs are the sources of ultra-high 

energy cosmic rays. 
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PACS numbers: 95.85.Ry, 14.60.Pq, 98.70.Rz, 98.70.Sa 



IceCube has become the most sensitive TeV-scale neutrino telescope that may reach the predicted neutrino flux 
from gamma-ray bursts (GRBs). The continued non-detection in its 22 (IC22, []J), 40 (IC40, 0) and 59 (IC59 Q) 
string configuration puts more and more stringent limits on GRB neutrino flux. The IC40 limit is comparable to 
the theoretical prediction Q, whereas the combined IC40 and IC59 limit is only 0.22 times the prediction Q. These 
limits start to put interesting constraints on the GRB neutrino models of Waxman and Bahcall [H, [1| and Guetta 
et al. (2004) Q, and challenges GRBs as the sources of ultra- high energy cosmic rays (UHECRs). However, in this 
note we show that the approach that IceCube [H-0| takes in theoretical prediction is somewhat different from that of 
Waxman and Bahcall and Guetta et al. (H-Qj leading to overestimate of GRB neutrino flux by a factor of ~ 5 for 
typical GRB parameters. This is because of ignoring the effects of the energy dependence of charged pion production 
and secondary pion/muon cooling on the normalization of neutrino flux. 

The approach taken by IceCube papers [IHH is presented in the appendix of the IC22 paper l]. In their approach, 
the muon neutrino flux from a GRB, T v , (with neutrino oscillation the electron, muon and tau neutrinos roughly 
share equal energy Q) is scaled to the proton flux in the GRB, T v , as 

o : T ^i^ = c 1 ) 

, (this comes from eq [A8] of [1]). Here = /-^(E = E b ) is the fraction of energy of protons with E b carried by 
charged pions, and E b is the energy of protons that interact with photons with spectral-break energy e b at A resonance 

a. 



E b = 1.3 x 10 16 r| 5 e-^ eV, (2) 



where T = 10 2 5 r 2 .5 is the bulk Lorentz factor of the GRB, and e b = le&,MeV MeV. The f„ :b value depends on the 
GRB properties, and varies from burst to burst (see eqs [A 7] and [A8] of [l|). The factor 1/8 is due to the facts 
that one half of the pj interactions produce charged pions, and that each generated neutrino is assumed to carry one 
fourth of the secondary pion energy. The proton flux can be normalized to gamma-ray flux by T p = (1/ f e )J-j where 
/ e is the ratio of accelerated proton to electron energy. In what follows we show that the approximation in eq ([l} is 
different from the models of Waxman and Bahcall [3,111 an d Guetta et al , and leads to overestimate of the neutrino 
flux. 

For a flat proton distribution with index p w 2 (E 2 dn p /dE cx E 2 ~ p ), and a typical GRB spectrum with Band- 
function parameters, a 7 = 1 and /3 7 = 2, f n (E) = is valid for protons with E > E b . However, for E < E b , 
fn oc E reduces with E decreasing because fewer target photons at high energy Q. Thus by using fir(E) — /„-,;, in 
energies E < E b eq. (UJ overestimates the neutrino flux. 

Moreover, eq. ([1} also i gno res the suppression of neutrino production at high energies due to the radiative cooling 
of secondary pions / muons [15J. The synchrotron cooling timescale is shorter than the secondary decay time at energies 
above the cooling energy [lUll [1] , 



E c>n/fl = 2 x 10 17 (e e /e B ) 1 / 2 rl 5 At_ 2 L- 1 / 2 x } n ^ eV. (3) 



S CiW/M = 2 x lU-(e e /e B r'-l'2.5^-2^52 1/2 x { lo (Sr±j 

Here L = 10 52 Ls2ergs _1 is the GRB (isotropic) luminosity, At — 10~ 2 At_2S is the GRB variability time, and e e and 
cb are the fractions of internal energy carried by postshock electrons and magnetic field, respectively. 
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Thus the neutrino production is mainly contributed by protons with E b < E < E Cl which is only a fraction of the 
total accelerated protons in energy. The distribution of the accelerated protons is expected to be a power law between 
the minimum and maximum energy. For mildly-relativistic GRB internal shocks, the minimum accelerated proton 
energy might be 



E mia ~ Tm p c 2 = 3 x 10 n r3ooeV. (4) 
The maximum proton energy is determined by the limit of synchrotron cooling for typical GRB parameters Q , 

£ max = 2.5 x 10 20 r^A^4 2 e y 4 eB 1 / 4 g- 1 / 2i -V4 eV; (g) 

where g > 1 accounts for the uncertainty in particle acceleration time. 

Since the neutrino production is mainly contributed by protons with Ef, < E < E c where f^(E) = f W) b, the muon 
neutrino flux T v is estimated to be 
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(6) 



where T p = /J' max E^g-dE and Eq Jl} have been used, and the last equation holds for p » 2. Given E max /E min ~ 10 9 
and E c ^/E b ~ 10 2 , the correction to Eq (QJ is a factor of F u /Fl c ~ 0.22. 

Below is more detailed calculation about TvlT^. For a flat proton distribution, p — 2, and a GRB spectrum with 
Band-function parameters, a 7 , /3 7 and e b , the source neutrino spectrum (before neutrino oscillation) is [H, for 
muon neutrinos from secondary pion decay, 



dnl 
de 



\{t) 



n x < 



for electron and muon neutrinos from secondary muon decay, 



e < e b 
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(8) 



There is no tau neutrino v T generated in pj interactions, and the small effect of kaon production on neutrino flux is 
neglected. Here the power law indices are a v = 3 — /3 7 and § v = 3 — a 7 , and the normalization is 
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(9) 



Note the factor 50 = \ x 0^/0.05 is resulted from the facts that (i) 1/2 of pq interactions produce charged pions; (ii) 
since the pion carries 0.2 of the proton energy, the muon neutrino number (after oscillation) per proton is 
and (iii) a single neutrino carries 0.2 X x = 0.05 of the proton energy, e b = 0.05E b . Similarly, we assume the other 
break energies in the neutrino spectrum, e Cy7T ~ 0.05.E C)7r , and e C)M ~ 0.05-E CjM . 

Considering neutrino oscillation, the muon neutrino spectrum (including and v>u) detected on the Earth is 
approximated as Q 
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Note here dn VT jde — since no tau neutrino produced. The muon neutrino flux is calculated as 

r 0.05E„ 



f 

JO. 



05B n 



: e(dn® /de)de 



jv*w E{dn p /dE)dE 



(10) 



(11) 



3 



For various values of GRB spectral parameters, a 7 , /3 7 and e^, we calculate the value of F v /F p with eqs. (|7] fTT|) . and 
show the T V IT„ value in Table HI We find that the approximation of eq ((T|) overestimates GRB neutrino flux by a 
few. Typically for GRBs with a 7 = 1, /3 7 = 2, and e£ bs = e b /(l + z) = 0.2 MeV, the correction factor is T v jT 1 ^ = 0.2, 
i.e., the neutrino flux is overestimated by a factor of 5 in [l|-|3j. 

Note the neutrino spectrum used here (eqs. [7] and [8]) is different from the one by IceCube papers [l|4j at high 
energy; they assume a simple steepening dn°^/de oc £ ~(f 3 »+ 2 ) at e > £ CijU (see eq. [A3] in IC22 paper which 
underestimate the neutrino emission from pion decay. Using their spectral shape, the neutrino flux is even smaller 
than using ours, thus the correction factor is even smaller. 



TABLE I: The value of T v /Tl c for various spectral parameter values. 



a 7 
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0.5 


1.5 
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1 1 


Pi 
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2 


2 


1.5 
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2 2 


ef s /MeV 


0.2 


0.2 


0.2 


0.2 


0.2 


0.05 2 




0.20 


0.12 


0.56 


0.25 


0.18 


0.14 0.30 



Note-The other parameters are L = 300, At = 10" 2 s, L = 10 52 ergs _1 , z = 1.5 and £? max = 10 21 eV. The redshift z only has 

effect on the GRB spectral break, = £t bs (l + z). 

In summary, for a flat energy distribution of accelerated protons in GRBs, p « 2, the predicted neutrino flux using 
eq (TJJ), so as the approach taken by [H-Q, is overestimated by a factor of J 7 * jT vl typically ~ 5, compared to the 
models by 0-Q . So the IC40 limit of neutrino flux Q actually should be above the prediction after correction, while 
the combined IC40 and IC59 limit, which is claimed to be 0.22 times of the prediction could be consistent with 
the correct prediction. 

Some comments should be made here. In this brief note we point out that when trying to test the models of 
with data, the IceCube papers [H-Q actually take an approach different from the former models in calculating 
the neutrino flux, leadin g to overestimate. However, the models of 0-0] consider only A resonance and neglect the 
multi-pion production P llOjl , the kaon production [ll[ and the possible secondary particle acceleration [l2| ■ All these 
facts may increase the neutrino flux and somewhat compensate the overestimate. It is also worth mentioning that if 
UHECRs are produced by the decay of neutrons that escape from the GRB outflow, it is straightforward to relate 
the neutrino flux with the observed UHECR flux, and avoid the uncertainties in parameters, e.g., /„■ and f e [l3j . If 
the assumption that only neutrons escape is true, the neutrino flux will also increase. 

Finally it should be stressed that there are large uncertainties in the GRB neutrino models, which may cause lower 
fn and hence lower neutrino flux, e.g., larger emission size or larger bulk Lorcntz factor. More observations by IceCube 
are required to push the sensitivity below the uncertainties [14(, and test the assumption of GRBs as the UHECR 
sources. 
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